Frustrated three-leg spin tubes: from spin 1/2 with chirality to spin 3/2 
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Motivated by the recent discovery of the spin tube [(CuCl2tachH)3Cl]Cl2 , we investigate the prop- 
erties of a frustrated three-leg spin tube with antiferromagnetic intra-ring and inter-ring couplings. 
We pay special attention to the evolution of the properties from weak to strong inter-ring coupling 
and show on the basis of extensive density matrix renormalization group and exact diagonalization 
calculations that the system undergoes a first-order phase transition between a dimerized gapped 
phase at weak coupling that can be described by the usual spin-chirality model and a gapless crit- 
ical phase at strong coupling that can be described by an effective spin-3/2 model. We also show 
that there is a magnetization plateau at 1/3 in the whole gapped phase and slightly beyond. The 
implications for [(CuCl2tachH)3Cl]Cl2 are discussed, with the conclusion that this system behaves 
essentially as a spin-3/2 chain. 

PACS numbers: 75.10.Dg,75.60.Ej 



I. INTRODUCTION 

During the last fifteen years, spin ladders have at- 
tracted a lot of attention both from theoretical and ex- 
perimental physicists (for an early review, see Ref. pj). 
On one hand, they provide a calculationally tractable 
system interpolating between one-dimensional and two- 
dimensional quantum spin physics. On the other hand, 
existing spin ladder compounds allow for a detailed com- 
parison of theoretical predictions and experimental re- 
sults. The interpolation between one and two dimensions 
is however not smooth: while half integer spin ladders 
with an even number of legs N show a gapped spectrum 
for spin excitations, ladders with odd N are believed to 
be gapless (like the spin-1/2 Heisenberg chain). 

For odd TV-ladders with half-integer spins, applying 
periodic boundary conditions in the transverse direc- 
tion of the ladder (forming thus a spin tube) yields an 
even more intriguing situation: In this case, for anti- 
ferromagnetic couplings, the ground state of an isolated 
ring is fourfold degenerate (twofold in spin- and twofold 
in chirality-space) and the low-energy physics of weakly 
coupled rings therefore involves both spin and chirality 
degrees of freedom. A number of theoretical studies have 
investigated effective low-energy Hamiltonians for weakly 
coupled rings by means of bosonizatior. 2,3 , density matrix 
renormalization group (DMRG) 4 '-, mean field theory^, 
and exact diagonalization (ED) 5 . All of them conclude 
that the ground state is dimerized, displaying a gap for 
both spin and chirality excitations. The details of the 
gaps depend however on the frustration of the inter-ring 
coupling 5 . 

Two experimental candidates for spin tubes are cur- 
rently available: The vanadium oxide Na2Va07i2iS, 
probably a nine-leg spin tube but with only three- 
fold rotational symmetry, and the three-leg compound 
[(CuCl2tachH) 3 Cl]Cl2i2iii. None of these examples is in 



the regime of weak coupling between the rings, intra- and 
inter-ring coupling constants being of the same order. It 
is therefore a very important first question to ask whether 
the low-energy physics of these realistic spin tubes may 
still be described by the above mentioned spin-chirality 
models. In both experimental realizations, neighboring 
rings are not coupled by one single antiferromagnetic 
bond, but by at least two competing bonds. This leads 
to additional frustration and raises the second question 
as to how much this frustration can change the physics. 

We address both questions by considering the exam- 
ple of a three-leg tube with antiferromagnetic intra-ring 
couplings and two frustrating antiferromagnetic inter- 
ring couplings which forms a minimal setup including 
chirality and additional frustration. This model repre- 
sents directly the compound [(CuCl2tachH)3Cl]Cl2, but 
the main results are expected to apply to other frustrated 
spin tubes as well. Regarding the properties, we will con- 
centrate on the presence or absence of a spin gap, on the 
nature of the low-lying excitations, and on the magneti- 
zation curve. 

The paper is organized as follows: In Sec. [H] we in- 
troduce the basic notations and describe several limiting 
cases of the three- leg spin tube. In Sec. lIIII we discuss the 
zero-field phase diagram from a numerical point of view 
and report on low-spin boundary excitations. In Sec. IIVI 
we use continuous unitary transformations to derive ef- 
fective low-energy Hamiltonians. Finally, we analyze the 
phase diagram as a function of the external magnetic 
field in Sec. Ivland discuss the experimental implications 
of our results in Sec. IVII 



II. THE MODEL 

The starting point is the Hamiltonian of the frustrated 
antiferromagnetic spin-1/2 Heisenberg model on a three- 
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FIG. 1: (Color on-line) Three-leg spin tube with ring cou- 
pling Ji (thick lines) and cross coupling J2 (thin lines). Both 
couplings are antiferromagnetic, leading to frustration both 
inside the triangles and between the triangles. 



leg spin tube: 



H 



Ji ai - 



J2 

ij 



(Oi+lj+l + Cj+lj-l) ■ 



(1) 

The index i runs along the tube, j = 1 ... 3 around the 
tube, and a is the spin operator for a spin 1/2. The tube 
can be viewed as series of triangles (see Fig. Each 
corner of such a triangle is coupled to two corners on the 
previous and the next triangle. In the absence of the J\ 
coupling, the Hamiltonian is bipartite (see Fig. m, and 
can be seen as a particular wrapping of a square lattice 
onto the tube. 



A. Lieb Schultz Mattis Theorem 

The Lieb Schultz Mattis (LSM) theorem^ is one of 
the few exact results on frustrated spin systems. Since 
its proof is quite general, it can be applied to a wide 
variety of systems with short-range interactionsi2ii4 ; in 
particular, to the spin tube Hamiltonian |T|. An explicit 
proof of the LSM theorem for this case can be found 
in Ref. [l^. The theorem states that half- integer spin 
chains have either a degenerate ground state or a gapless 
spectrum due to a zero-energy mode with momentum tt 
relative to the ground state. It is interesting to note that 
both situations are realized in the three-leg spin tube 
system. We first discuss them in two limiting cases. 



B. Weakly coupled triangles 

In the limit J\ ^> J 2, the system consists of weakly 
interacting triangles. Each isolated triangle has a four- 
fold degenerate spin- 1/2 ground state that can be labeled 
by its chirality, r z , and its magnetic moment along the z 
axis, S z . The projection of the full Hamiltonian Q onto 



the low-energy space spanned by the spin-1/2 states leads 
to lowest order in J 2 to an effective spin-chirality model, 

H1/2 = Ps=i/2HPs=i/2 



2J2 
3 



■E i S i -S i+ i(l + a[r+.r i ; 1 +h.c.]) , 

(2) 

where a = 2 for the frustrated tube of Eq. QJ. Here, Sj 
is a spin-1/2 operator describing the total spin on trian- 
gle i, Ti is a pseudo-spin- 1/2 operator acting in chirality 
space, and Ps=i/2 is the projector in the subspace where 
each triangle is in a S = 1/2 state. It can be written as 



the product of the local projectors, Ps=i/2 
with: 



= n,p, 



p, 



S=l/2 



15 - 4S| 
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S=l/2> 



(3) 



The reduced spin-chirality Hamiltonian J5J has been 
studied in Refs. 00. For a = 2, it displays a dimer- 
ized ground state with a twofold degeneracy and gaps 
for both spin and chirality excitations. 



C. Strongly coupled triangles 

In the limit J\ = 0, the lattice is effectively bipartite, 
and we expect ferromagnetic correlations between spins 
on the same ring because they belong to the same sublat- 
tice. The total spin on each ring should therefore be close 
to its maximum value of 3/2, and the physics for Ji = 
should be similar to the case of a ferromagnetic coupling, 
Ji < 0, within each triangle. The effective Hamiltonian 
in the strong coupling case —J\ 3> J2 reads: 



H 



3/2 



(4) 



where St is now a spin-3/2 operator. The coupling 
constants are given by K\ — 2J 2 /3 + jf/(18|Ji|) and 
A'2 = — Jf/(54| Ji|). Note that the nature of the second 
order corrections to K\ and K2 does not seem to frustrate 
the effective spin-3/2 chain. It is therefore reasonable to 
argue that the spin-3/2 state is stabilized at least up to 
Ji = even though this value lies beyond the pertur- 
bative limit. This argument can be made more quanti- 
tative by performing range-two CORE-calculationsiSiii 
yielding effective coupling constants K\ « 0.78 and 
K2 ~ —0.04 for J 1 — (the bicubic K3 coefficient is 4 
times smaller than 1^2 )• The ratio of the couplings does 
not change drastically even by increasing J1/J2 up to 1. 
The physics of the effective Hamiltonian J2J is therefore 
dominated by the Heisenberg (Ki) term, which is known 
to belong to the universality class of a Luttinger liquidiS, 
and therefore has a gapless spectrum. We believe that 
small biquadratic or bicubic interactions will not imme- 
diately open a gapi^ and we expect the low energy theory 
to remain stable. Note that these gapless excitations are 
in sharp contrast to the gapped spectrum of the spin- 
chirality Hamiltonian J2J). 
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Based on the preliminary arguments of subsections lll Bl 
and 111 Cl we expect the frustrated spin-tube Hamiltonian 
l|T| to undergo at least one quantum phase transition 
as we move from weakly coupled ( J\ 3> J2) to strongly 
coupled {J 1 <C J2) triangles. 



III. NUMERICAL RESULTS 

In the following, we investigate the properties of this 
phase transition numerically. We have performed ex- 
tensive density matrix renormalization group (DMRG) 
calculations 20 on tubes with open boundaries and L < 
100 and exact diagonalization (ED) on small systems up 
to L = 12 using periodic boundary conditions (PBC). 

Two independent quantities prove clearly that there is 
a single first-order bulk phase transition. They are dis- 
cussed in Secs. lIII Al and lTTl Bl respectivelv. Furthermore, 
on open boundary systems, there is a second "transition" 
at which the state of rings at the boundary jump from 
predominantly S = 1/2 to predominantly S = 3/2. This 
phenomenon will be discussed in Sec. lIIICl 



A. Local spin on a ring 

The first-order character is best seen in the ground 
state expectation value (Ps =1 / 2 ) °f the previously defined 
projector ©, which is a purely local quantity defined on 
a single triangle We have computed this expectation 
value as a function of J2. The results are shown in Fig. [3 
which displays the projector in the middle of the tube at 
i = L/2. For weak coupling, J2/J1 <C 1, the expectation 
value is close to one, for strong coupling, J2/J1 ~> 1, it 
is close to zero. The expectation value clearly jumps at 
a value J 2c = (1.219 ± 0.003) J\ indicating a first-order 
phase transition which is magnified in the inset of Fig. [21 

Several comments are in order: Fig. [21 displays results 
obtained both with DMRG and ED. To check the im- 
portance of boundary effects (see Sec. lIIICl) . DMRG was 
performed with two types of boundary conditions: the 
standard Open Boundary Conditions (OBC), and the 
Ferromagnetic Boundary Conditions (FBC), where we 
put a ferromagnetic coupling J\ = —10 on the two tri- 
angles at the boundary of the tube, thereby strongly fa- 
voring the 5 = 3/2 state. For the ED results, periodic 
boundary conditions were applied. The curves for the 
three kinds of boundary conditions differ only near the 
first-order transition, and the differences decrease with 
the increasing length of the tube. In the ED, we addi- 
tionally observe that for system sizes L = 4p + 2 the 
ground state changes the momentum sector at the tran- 
sition from k = for weak coupling to k — tt for strong 
coupling. This is another fingerprint of a first-order tran- 
sition. 
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FIG. 2: (Color on-line) Expectation value of the projector 
m t ne bulk as a function of J2/J1 for different lengths 
and boundary conditions (see text). The jump at J%jJ\ fa 
1.22 is a clear indication of a first-order transition. 
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FIG. 3: (Color on-line) Dimerization at the center of the tube 
Od as a function of J2/J1 for L = 50 (black squares) and 
L — 100 (red circles). The vertical line represents the dis- 
appearance of the true order parameter: The dimerization 
above this line vanishes after proper finite-size scaling. 



B. Dimerization 

The first-order character of the phase transition is also 
clearly apparent in inter-ring correlations. As an exam- 
ple, we investigate the local dimerization defined as 



Di — (— I) 1 (Sj.Sj+i — Si + i.Si+2) . 

This dimerization can be viewed as the order parameter 
of the symmetry-broken weak coupling phase. It is also 
present in the spin-chirality model discussed in Sec. Ill Bl 
for which the dimerization opens a spin gap^*^. Since 
we work with open boundary conditions, the quantity Di 
varies with the ring position i. Fig.[31shows the order pa- 
rameter Od{L) = (-D.L/2) (the dimerization in the middle 




FIG. 4: (Color on-line) Spatial dependence of the projec- 
tor {-Pg = i/2) (upper panel) and nearest-neighbor correlation 
function (Si ■ Si+i) (lower panel) for L — 50 and different 
values of the inter-ring coupling J2 . 

of the tube) for different system sizes L. Two sharp tran- 
sitions are observed, the first one around J2 = J2 C and 
the second one at a higher value Jij J\ ~ 1.47. How- 
ever, finite-size scaling shows that the order parameter 
Od{L) remains finite only for J2 < J2 C in the thermo- 
dynamic limit. The first transition, the disappearance of 
the dimerization, corresponds to the phase transition to 
a 5 = 3/2 phase on the rings. The second transition is a 
boundary effect discussed in the next section. 



C. Boundary excitations 

Boundary spin- 1/2 degrees of freedom in open spin 
5 = 3/2 chains have been predicted theoretically in 
Ref. H, and were later confirmed in DMRG studies of 
unfrustrated^ and frustrated^ spin 5 = 3/2 chains. 

In our model we find these edge states as well for suffi- 
ciently large JiJ J\. However approaching the first-order 
transition coming from J 2 > J20 there is a second class 
of edge states, which are related to a kind of nucleation 
of the dimerized 5=1/2 phase at the boundaries. These 
edge states arc different than those discussed above, as 
they now originate from the 5 = 1/2 subspace of a ring, 
and therefore also include a chirality degree of freedom. 

This is most clearly seen in the spatial dependence of 
the correlation functions considered in the two subsec- 
tions above. The upper panel of Fig. 0] shows the space 
dependence of the projector (P§ =1 / 2 ) and the lower panel 
of the nearest-neighbor spin correlation (Sj • Sj-|_i) for dif- 
ferent values of J<zj J\. The ring at the boundary stays 
in a predominant 5=1/2 state even for J2 > J 2c , where 
the bulk transition has already occurred. The boundary 
spin- 1/2 states disappear only beyond J 2 ~ 1.5. The 
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FIG. 5: (Color on-line) The projector (P| =1 / a ) in the bulk 
and at the boundary of the tube for L — 100 jumps for differ- 
ent critical couplings indicating the existence of a boundary 
excitation. 

nearest-neighbor spin-spin correlation function follows a 
similar behavior: the correlations at the boundary are 
close to those obtained in the bulk 5 = 1/2 phase, and 
then gradually approach the values of the bulk 5 = 3/2 
phase. 

The different behavior of bulk and edge rings are nicely 
illustrated in Fig. in which the values of {P(s=i/2)) i n 
the bulk and at the boundary of the tube arc plotted as 
a function of the coupling ratio J2/ J±. The expectation 
value of the projector drops at J2/J1 ~ 1.22, whereas 
the boundary ring remains in a 5 = 1/2 state up to 
Ja/Jiwl.5. 

IV. EFFECTIVE HAMILTONIAN APPROACH 

Since the numerical results strongly indicate that the 
transition is first order, it should be possible to describe 
the transition as occuring directly between an effective 
spin-chirality model and a spin 3/2 model. Indeed, we 
have previously discussed that the ground state of the 
model Hamiltonian lies in the spin-1/2 sector for 
J\ » J2 and in the spin-3/2 sector for J 2 <C J\. The 
problem is then reduced to diagonalizing an effective 
Hamiltonian @ or jl| on a truncated Hilbert space (for 
the spin-1/2 case see Refs. QUI). However, in the in- 
termediate regime, J\ ~ J2, it is no longer possible to 
construct effective Hamiltonians using perturbation the- 
ory in the small parameter J2/ \ J\\ and an alternative 
construction has to be found. In this context, it is inter- 
esting to observe that the above Hamiltonian H can be 
extended to 

H e = H + J' 2 ^ (Tij ■ <T i+ i j. 

ij 

For J2 — J2, the extended Hamiltonian H e has the spe- 
cial property to conserve total spin and chirality of each 
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FIG. 6: (Color on-line) Comparison of ground state energies 
in the spin-1/2 and spin-3/2 sectors as a function of J2/J1 
with 15 couplings kept in the flow. Finite-size effects only 
appear for short tubes of length N = 6,8. The first-order 
phase transition is located near J2 = 1.4 Ji. The vertical line 
represents the position of the phase transition according to 
the numerical results. The discrepancy between the two is 
due to the truncation in the space of possible couplings. 



sired properties is obtained 
dH(£) 



df 



Hi 



00 



H, 



cir- 



Following Ref. |26j], we choose the commutator 

n{i) = [H c {i),H{i)] 



(5) 



(G) 



as an appropriate generator of the transformation. Here, 
H c (£) denotes the conserving part of the Hamiltonian 
which does not change the number of spin-1/2 triangles. 
The above defined flow generates an infinite set of cou- 
plings which has to be suitably truncated. To illustrate 
the behavior of the flow qualitatively, we first use a trun- 
cated set of only three couplings. Later on, we will use 
a set of sixteen couplings which allows us to arrive at 
semi-quantitative results. 

For the qualitative analysis, we decompose the Hamil- 
tonian 

H(£) = Jxi^hx + J %c {l)h %c + J 2 , n (£)h2, n , 
into the coupling 

hi = <7jj • (Tij + l, 



triangle separatel}^. Note that this property is inde- 
pendent of the relative strength of J\ and J2 and holds 
also for Ji ~ J 2 . The Hamiltonian H e is naturally block- 
diagonal when J2 — J' 2 , each block being characterized 
by a set of quantum numbers for spin and chirality on 
each triangle. The ground state is found either in the 
block where all triangles have spin-1/2 or where all have 
spin-3/2. Both cases are obviously connected by a first 
order phase transition. 

For our Hamiltonian H we have J2 ^ J'2 = and the 
large set of conserved quantum numbers is lost. Even the 
total number of spin-1/2 triangles is no longer conserved. 
To restore this latter conservation law, we propose to ap- 
ply a continuous unitary transformation to the Hamilto- 
nian H which eliminates all couplings that change the 
total number of spin-1/2 triangles. The obtained effec- 
tive Hamiltonian H c g is again block-diagonal and allows 
us to calculate the ground state energies of the spin-1/2 
and spin-3/2 sector independently. We may then infer 
the critical J2 C for the first order phase transition from 
the comparison of those energies. The proposed unitary 
transformation can be carried out exactly in the limit 
J 2 — J' 2 — > 0. In our case, J 2 = 0, the construction of the 
effective Hamiltonian H c g can be viewed as a perturba- 
tive expansion in the parameter J2 — J 2 . 



A powerful tool for the derivation of effective Hamil- 
tonians was introduced by Wegner— , the so-called flow 
equations which were first applied to spin Hamiltoni- 
ans by UhrigSi. Flow equations apply infinitesimal uni- 
tary transformations (parametrized by £) to the original 
Hamiltonian until an effective Hamiltonian with the de- 



the coupling ft,2 jC which is the part of 

ll2 = " (O'i+lj + l + O'i+lj-l) 

ij 

that conserves the number of spin-1/2 triangles and the 
coupling h2.n = h 2 — /i2,c that changes the number of 
spin-1/2 triangles. The two couplings ft.2,c and ft.2,n are 
readily obtained by applying the projectors (0 to the 
terms of /12. 

Inserting the Hamiltonian H (£) and its conserving part 
H c (£) = Jt^hx + J 2 ,c(^)/i2, c into Eqs. © and ®, we 
obtain (after the truncation) three coupled differential 
equations for the coupling constants 



dJi 
~d£ 



32Ji -2I-J2,, 



dJ- 



2.c 



df 



K in 
13yJi + 19— J 2 , c ) (J 2 ,„) 2 



dj2 > n = f -72J 2 



( // : 



1 . t 96JiJa,c-67jjJ£ c ) J a ,n, (7) 



with the initial conditions 

Ji(0) = Ji,J 2 , c (0) = J a ,„(0) = J 2 . 

The analysis of the last differential equation shows that 
J2.n scales to zero during the flow, leading to an effective 
Hamiltonian H c fi which obeys the desired conservation 
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law. The qualitative inspection of the first and second 
equation shows that for J 2 < 1.5Ji the spin-1/2 phase is 
stabilized (Ji(oo) > Ji(0), J2, c (oo) < J2, c (0)), whereas 
for Ji > 1.5Ji the spin-3/2 phase is stabilized in the 
same way. For very large initial Ji the effective Ji(oo) 
becomes even ferromagnetic Ji(oo) < 0, so that the one- 
site coupling hi itself favors the spin-3/2 state. This 
analysis is confirmed by the numerical integration of Eqs. 
0. 

The accuracy of the above reasoning can be made semi- 
quantitative, when more couplings are included in the 
flow. For this purpose we evolve the one-triangle cou- 
pling hi and all the fifteen two-triangle nearest-neighbor 
couplings which are allowed by symmetry under the flow. 
We evaluate the ground state energies in the spin-1/2 
and spin-3/2 sectors of H c g by exact diagonalization of 
small chains corresponding to tubes up to length L = 14. 
(Note that the dimension of the effective Hilbert space 
is reduced in comparison to the full Hilbert space by a 
factor of 2 L ). The results of the numerical calculation 
are shown in Fig[5] The phase transition between a spin- 
1/2-chirality chain and a spin-3/2 chain appears around 
J 2 ~ 1-4 J\. This value lies slightly above the precise J 2c 
obtained by DMRG. As can be seen from Fig. finite- 
size effects do not explain this quantitative discrepancy. 
It is in fact due to the truncation of the flow equations 
to nearest-neighbor couplings. It can be checked numer- 
ically that the effect of including next-nearest neighbor 
two-triangle couplings in the flow is negligible. To reach a 
more quantitative result, it is necessary to include three- 
triangle couplings. 

In summary, a first-order transition between a spin- 
chirality model and a spin-3/2 model is also found in 
this approach, and the agreement with the critical ratio 
J2/J1 found in DMRG is very satisfactory considering 
the truncation of the flow. Therefore we believe that the 
case for a first-order transition is very strong. 



V. PHASE DIAGRAM IN A FIELD 

Finally, to make contact with the experiments per- 
formed on [(CuCl2tachH)3Cl]Cl2, we investigate the 
properties of the model in the presence of an external 
magnetic field. 

For J2/J1 = 0, the triangles are decoupled and we ex- 
pect a magnetization jump from m = to to = 1/3 at 
h = + , where all the triangles are in a spin-1/2 state 
with S z = +1/2. This to — 1/3 plateau should survive 
in a finite parameter range. In terms of the spin-chirality 
model of Ref. [j| , such a plateau is quite unusual. Indeed, 
when the spin degrees of freedom align ferromagnetically 
in the field, the effective model for the chirality becomes 
an XY model, whose spectrum is gapless. So, accord- 
ing to this argument, this plateau phase is expected to 
have a gap to spin excitations - as do all plateau phases 
- and to simultaneously possess gapless chirality excita- 
tions. Note that higher-order terms in the effective model 
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FIG. 7: (Color on-line) Magnetization curves for a spin tube 
with L = 36 and three different values of Jij J\. In the right 
panel, we also show the magnetization curve of a S = 3/2 
chain for comparison. 

might open a small gap in the chirality excitations as well, 
but being a higher order effect, this gap (if any) should 
be small compared to the spin gap, and should thus lead 
to observable properties, such as a specific heat linear in 
T below the temperature scale set by the width of the 
plateau. 

Typical magnetization curves are shown in Fig. [7| for 
several values of J2/J1. As long as the ratio J2/J1 is not 
too large, a plateau at m = 1/3 is clearly present. As 
can be seen in Fig. [HI this plateau phase survives every- 
where in the region in which the system is effectively a 
spin-1/2 chain with chirality, but interestingly, it extends 
beyond that point up to J 2 /Ji ~ 1.6. It is plausible that 
the region of validity of the effective spin-l/2-chirality 
model extends to higher values of J2/ Ji when a magnetic 
field is applied to the system. Indeed, an abrupt jump 
in Pg^i/2 i s stm present for slightly polarized samples 
and the corresponding critical value of the ratio J2/J1 
increases slightly with magnetization. By the time the 
m = 1/3 plateau is reached, the suppression of Pg—i/% is 
very smooth, however, and it is not possible to decide on 
the basis of the available data whether the first-order line 
extends up to the right boundary of the plateau phase, 
or whether it ends at a critical point on the way. In any 
case, for large J2/J1, the magnetization curve becomes 
smooth, as it should for a spin-3/2 chain. 



VI. [(CuCl2tachH) 3 Cl]Cl 2 

The Hamiltonian of Eq. Q is believed to be realized in 
the compound [(CuCl 2 tachH) 3 Cl]Cl 2 , with J 2 /Ji » 2.16, 
and the theoretical investigation of that model reported 
in Ref. \\\\ came to the somewhat surprising conclusion 
that the system has a gap to all (spin and singlet) excita- 
tions. This is at odds with the LSM theorem (see section 
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FIG. 8: Phase diagram as a function of J2/J1 and of the 
magnetic field H. Symbols are numerical data and solid lines 
are guides to the eye for the boundaries of the phase transi- 
tions. The gray area is a gapless phase, and the white area 
are plateau phases with a spin gap. Dotted line: experimental 
value of J2/J1 for the compound of Ref.[Tl| 
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1/3 plateau, there should be no plateau at 1/3. In fact, 
according to the discussion of the previous section, the 
absence of a 1/3 plateau in the magnetization curve of 
[(CuCl2tachH)3Cl]Cl2 reported in Ref. ^lj shows unam- 
biguously that this compound lies in the effective S = 3/2 
strong coupling phase. 

To be more specific, the claim regarding the absence of 
a spin gap is clearly supported by our calculations of the 
gap to S — 2 excitations (since we are working with open 
boundary conditions, the gap to S = 1 excitations is very 
small and corresponds to boundary excitations 22 23 ). In- 
deed, the finite-size scaling of the gap is consistent with a 
vanishing value in the thermodynamic limit for all ratios 
J2/J1 larger than 1.22 (see Fig.|5J). 

Regarding the magnetization curve, our results agree 
with those of Ref. 11], which were obtained on smaller 
systems: The curve is smooth with no trace of any ir- 
regularity close to 1/3-magnetization. Fig. [7| shows the 
magnetization curve for the spin tube for various value 
of J2 and for the spin-3/2 chain for L = 36. The cou- 
pling of the spin-3/2 chain has been chosen so that the 
saturation field is the same as the saturation field for the 
experimental value of the couplings (J2/J1 = 2.16). The 
two magnetization curves are very similar giving further 
support to our claim that [(CuCl2tachH) 3 Cl]Cl2 should 
be regarded as a spin-3/2 chain. It would be desirable 
to have additional low-temperature susceptibility data 
in order to clarify whether the compound really shows 
a spin gap. On the basis of the available susceptibility 
and magnetization data, gapless behavior is clearly not 
excluded. 



VII. CONCLUSIONS 



°0 0.005 0.01 0.015 0.02 0.025 0.03 
1/L 

FIG. 9: (Color on-line) Gap to the S=2 sector as a function 
of 1/L for several values of J2/J1 



III Aj) . On the basis of our DMRG results, we do not reach 
the same conclusion. If J2/J1 is small enough, the sys- 
tem should be in a dimerized phase with a spin gap and 
a twofold degenerate ground state. Hence a low-lying 
singlet excitation should be present in finite systems and 
should collapse onto the ground state in the thermody- 
namic limit. If J2/ J\ is large enough, the system should 
essentially behave as a spin-3/2 chain and be gapless with 
low- lying spinon excitations. The phase transition be- 
tween these phases takes place for J2/J1 « 1.22. So, 
with a ratio J2/J1 ~ 2.16, as deduced from the temper- 
ature dependence of the susceptibility in Ref. we 
predict that [(CuCl 2 tachH) 3 Cl]Cl2 should have a gap- 
less spectrum, and, since this ratio is also larger than 
the critical value that marks the disappearance of the 



We have investigated three-leg spin tubes consisting of 
rings that are coupled by competing antiferromagnetic 
bonds. We have shown that this frustration can drive 
the system away from the typically used effective spin- 
1/2-chirality model of regular spin tubes. It becomes an 
effective spin-3/2 chain even if the inter-ring couplings 
are antiferromagnetic. In the specific case of the model 
relevant to [(CuCl2tachH)3Cl]Cl2 in which each spin of 
a ring is coupled to two spins of the neighboring rings, 
the transition between these phases as a function of the 
inter-ring coupling has been shown to be first order. This 
should be contrasted to the case of the unfrustrated spin 
tube, in which case the transition is expected to be sec- 
ond order and to take place when the inter-ring coupling 
changes sign. These results lead to a new interpretation 
of the properties of [(CuCl2tachH)3Cl]Cl2, and to specific 
predictions regarding the occurrence of a plateau at 1/3 
in this geometry. It is our hope that the present work 
will motivate further investigations of the excitations of 
[(CuCl2tachH)3Cl]Cl2 to check our predictions, and more 
generally of the (so far) small but fascinating family of 
spin tubes. 



8 



Note added 

After submission of this work, K. Okunishi et al. in- 
formed us that they had submitted a work reaching the 
same conclusion concerning the existence of a first order 
phase transition^. 



by the Swiss National Fund and by MaNEP. Calculations 
have been partly performed on the Pleiades cluster of 
EPFL. 
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